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A MARGINAL DIMENSION OF A WEAKLY DILUTED 
QUENCHED m- VECTOR MODEL 
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We calculate a marginal order parameter dimension m c which in a weakly diluted quenched m- 
vector model controls the crossover from a universality class of a "pure" model (m > m c ) to a new 
universality class (m < m c ). Exploiting the Harris criterion and the field-theoretical renormalization 
group approach allows us to obtain m c as a five-loop e-expansion as well as a six-loop pseudo-e 
expansion. In order to estimate the numerical value of m c we process the series by precisely adjusted 
Pade-Borel-Leroy resummation procedures. Our final result m c = 1.912 ± 0.004 < 2 stems from 
the longer and more reliable pseudo-£ expansion, suggesting that a weak quenched disorder does 
not change the values of xy-model critical exponents as it follows from the experiments on critical 
properties of He 4 in porous media. 
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I. INTRODUCTION 

Along with the space dimension d, the order parame- 
ter dimension m is relevant for the universal properties 
of a model at criticality Q. In the presence of disor- 
der the value of m determines also whether a disordered 
model possesses novel universal properties in compari- 
son with a pure model. A typical example of this feature 
is presented by a critical behaviour of a weakly diluted 
quenched m-vector model [Q. The model is defined by a 
Hamiltonian 



H — — — ^ J(|Rj — Rj|)S'R i S'R j CR i CR j , 



(1) 



ho 



ity: the corresponding critical exponents remain negative 
a = -0.011 ± 0.004, a = -0.122 ± 0.009 §]. Therefore 
one can expect that only weakly diluted quenched Ising 
model belongs to a new universality class. 

Indeed, experimental studies confirm the theoretical 
prediction. The bulk of evidence collected in a recent re- 
view H demonstrate a novel critical behaviour of the 
magnetic systems described by model (Q) at m = 1. The 
experimental value of the heat capacity critical expo- 
nent at the A-transition in He-4 0] corroborated that 
the system belongs to the xy-model universality class 
with no divergency of the heat capacity at the phase 
transition point. Subsequently, experiments on the criti- 
cal behaviour of He-4 in porous media Q confirmed the 
irrelevancy of a weak quenched disorder in this case. 



where Ri span over the sites of a simple cubic lattice 
and Sr.; denote the m-component spins interacting via 
a translationally invariant short-range isotropic interac- 
tion J(|Ri — Rj|). The weak disorder is introduced by 
stochastic uncorrelated occupation numbers cr ; equal to 
1 in the case when a site is occupied by a magnetic atom 
and otherwise (see Fig. 1). The concentration of oc- 
cupied sites is considered to be above the percolation 
threshold. The quenched disorder implies that vacancies 
CR i = are fixed and require configurational averaging 
of observables M. 

The crucial dependence of the universality class of 
model (Q) on the order parameter dimension m can be 
established by the Harris criterion It states that a 
disorder changes the universal critical properties of a 
"pure" model only if heat capacity critical exponent of a 
pure model is positive. Within the hierarchy of the phys- 
ical realizations of the m-vector model only the Ising 
model m = 1 is characterised by a = 0.109 ± 0.004 > 0, 
while xy- (m = 2) and Heisenberg (m = 3) models' 
heat capacity does not exhibit divergency at critical- 




Fig. 1. The weakly diluted quenched m-vector model de- 
scribes a system of randomly distributed m-component vec- 
tors which are fixed on sites of a three-dimensional cubic lat- 
tice and interact via a short-range translationally-invariant 
force. 
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A natural question arises: can a marginal value of m 
be obtained such that for m > m c the critical exponents 
of a weakly diluted quenched m-vector model coincide 
with the corresponding values of a "pure" model, but for 
m < m c they split to new values. Note, that while the 
values of the critical exponents both for pure || and di- 
luted m-vector || models are calculated precisely, the 
crossover between two universality classes has not been 
the subject of extended theoretical calculations. One can 
mention only the results for m c on the basis of two- JlO| 
and three-loop renormalization group studies [pd| of a 
weakly diluted <i-dimensional m-vector model. The ob- 
tained values m c = 2.01 and m c = 2.12 contradict both 
experimental data ||] and the Harris criterion Q| . In con- 
trast, an alternative estimate on the base of a refined by 
conformal mapping resummation of the six-loop renor- 
malization group functions of the m-vector model ex- 
plicitly yielded m c = 1.942 ±0.026 flf]. In this study we 
want to solve the inconsistency of the theoretical results 
10 1^] and to perform an independent determination of 
m c . 

Relying on the Harris criterion one can determine 
m c from the requirement of vanishing the heat capac- 
ity critical exponent a of the m-vector model. We will 
make use of this condition appealing to a well-established 
field-theoretical renormalization group (RG) approach 
(lj which allows perturbative calculating of the marginal 
order parameter dimension m c . In the next Section we 
obtain m c as e- and pseudo-e expansions which follow 
from the alternative minimal subtraction and massive 
schemes of field-theoretical renormalization group ap- 
proach. In order to cope with the divergencies of the 
expansions and to obtain reliable numerical values based 
on them we process the expansions by appropriate re- 
summation procedures. The resummation procedures are 
also analyzed and adjusted in the next Section. The out- 
comes of the study are presented in the concluding Sec- 
tion 



III 



II. THE EXPANSIONS AND THE 
RESUMMATION 

As mentioned above, the marginal dimension m c of 
a weakly diluted m-vector model can be reconstituted 
from the critical properties of a "pure" m-vector model. 
Considering the heat capacity critical exponent a of a 
"pure" model as a function of the order parameter com- 
ponent number m, one can write the master equation for 
determining m c as follows: 



a(m c ) = 0. 



(2) 



The treatment of Eq. (g) by means of the field-theoretical 
RG approach can be performed in various schemes. Two 
of them are most widely used in the theory of critical phe- 
nomena. The dimensional regularisation with the mini- 
mal subtraction (l3j | allows obtaining quantities of in- 
terest by familiar e-expansions jl4| with e = 4 — d. As a 
starting point to obtain the e-expansion for m c serves the 



e-expansion for the critical exponent a of the m-vector 
model, which is known in the five-loop the approxima- 
tion p5[ . By keeping the coefficients of the expansions 
as functions of m and reexpanding equation (^|) in e we 
obtain m r in the form: 



4 - 4e + 4.707199e 2 - 8.727517e 3 



-20.878373e 4 



(3) 



Formally, the numerical value of m c at d = 3 can be cal- 
culated from the expansion (^) by the substitution e = 1. 
The renormalization conditions of the RG massive 



scheme |16 provide another possibility to obtain the crit- 
ical properties directly at d = 3. The traditional calcula- 
tion of the critical exponents in the massive field theoret- 
ical RG scheme implies a numerical analysis of the RG 
functions. However, the most accurate estimates of the 
critical exponents of a three-dimensional m-vector model 
P.[l7| are based on a pseudo-e expansion technique Jl8[ . 
This technique avoids intrinsic errors accumulation typi- 
cal for the numerical processing of the massive RG func- 
tions and results in a self-consistent collection of con- 
tributions from the different steps of calculations fll7j . 
On the other hand, introducting the auxiliary pseudo- 
e parameter r which mimics the role of the ordinary e 
parameter of the minimal subtraction RG scheme allows 
to analyse the series by methods well-established for the 
e-expansion. 

In order to obtain a pseudo-e expansion for m c we 
start from the expansions of RG functions of m-vector 
model which within the massive scheme have been cal- 
culated in the six-loop approximation ||19||. Introducing 
into functions |l9| the parameter r (see ]17| for details) 
we obtain the pseudo-e expansion for m c at d = 3 as 
follows: 



4 - 8/3t + 0.766489t 2 - 0.293632r 3 



-0.193141t 4 - 0.192714t 5 . 



(4) 



Again, the resulting numerical value of m c can be ob- 
tained from expansion (^) by the final substitution r = 1. 

The explicit form of expansions (||) and ([!]) is suffi- 
cient in principle to estimate the numerical value of m c . 
However, the series for RG functions are known to be of 
asymptotic nature pp|-p2]| and must be resummed before 
the final substitutions e = 1 (r = 1). An explicit form of 
the asymptotic of series (|)-(§ has not been obtained, 
contrary to e-expansions for the m-vector model critical 
exponents. Consequently, the applicability of resumma- 
tion procedures to the series of m c is only conjectured. 

A resummation procedure, which in different modi- 
fications is commonly used in the studies of asymptotic 
series, is the integral Borel transformation 23 1. However, 
this technique implies explicit knowledge of the general 
term of a series and thus cannot be applied here, where 
only truncated sums of the series are known. On the 
other hand neither can we use the resummation proce- 
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dures based on the conformal mapping technique |17j ] 
since even estimates on large-order order behaviour of 
expansions terms are unavailable. In the case of the se- 
ries (H)-(Q) we are restricted to the simplest procedures 
which do not imply such estimates. 

Let us start the analysis by representing series 
@ by means of Pade approximants [M/N] (x) — 
J2iLo a i x% / Ylj=o bjX-i in the variable x = e |Q. The 
result for m c is shown in the form of a Pade table (|5|), 
where the number of the row, M, and that of the column, 
N, correspond to the order of the numerator and that of 
the denominator of the Pade approximant [M/N] respec- 
tively. Subsequently, o denotes approximants which can 
not be constructed within the considered approximation, 
while small numbers in (^) indicate that approximants 
have poles (at e = 6.52 and e = 1.11 for [0/2] and [0/4] 
respectively) and thus are unreliable. 



six- loop approximants [3/2] and [2/3] as well as the five- 
loop approximant [2/2] yield practically the same value 
of m c within the fourth digit leading to the estimate 
m c = 1.9085. This number is not yet considered as the 
most accurate estimate of m c that one can obtain from 
expansion (Q). 

By assuming the factorial divergence of the coefficients 
of expansions (||), ([I]), one can apply to their analysis a 
refined Pade-Borel-Leroy resummation procedure which 
has been successfully used in various tasks of theory of 
critical phenomena pf|. The procedure is based on the 
integral Borel transformation |23| , however it uses as an 
intermediate step an extrapolation by means of a Pade- 
approximant |24| . More precisely, the procedure is de- 
fined by the following algorithm: 

• starting from the initial sum S of L terms one con- 
structs its Borel-Leroy image 



4 2 2 1939 1.5086 6 1528 

2.1624 2.0316 1.9365 o 

4.7072 1.6493 1.9208 o o 

-4.0203 2.1344 o o o 

16.858 oooo 



(5) 



The first column of table (|5|) represents a straightforward 
summation of series (|3|) terms and obviously shows the 
divergence of results with the increase of the approxi- 
mation order. Conversely, the convergence of numbers is 
observed along the main and next to main diagonals of 
the table. For instance, the expected inequality m c < 2 
can be obtained already within the fourth and the fifth 
order of the e-expansion (|J) . 

The results of Pade-analysis of the pseudo-e expansion 
(Q) are presented in table in the same notations as in 
table (|). 



4 2.4 2.0839 1.9669 1.9398 1.9106 

1.3333 1.9287 1 8799 1.9311 2 2425 o 

2.0998 1.8875 1.9084 1.9085 o o 

1.8062 1.9227 1.9085 o o o 

1.9993 1.9029 oooo 
1.8066 ooooo 



(0) 



Again, small numbers correspond to unreliable approx- 
imants with the poles at r = 7.29 and r = 0.907 for 
[1/2] and [1/4] respectively. The distinguishing property 
of table (H) is a convergence of the results on the basis 
of a mere summation of the pseudo-e expansion terms 
(Q) (the first column of table (^)). Contrary to the e- 
expansion analysis (||), the inequality m c < 2 is observed 
already in a four-loop approximation and remains valid 
in the five- and the six-loop approximations. However, 
the best convergence is noticed for the approximants par- 
allel to the main diagonal of table (|^). For instance, the 



s( x ) = J2 a ^ sB ( xt ) = 



4=0 



aijxt) 1 
(i + b)\ 



(7) 



where b is an arbitrary non-negative number; 

• subsequently one extrapolates the Borel-Leroy im- 
age (fil) by a rational Pade approximant 



S B (xt) 



[M/N] (art); 



• the resummed function S Res is finally obtained in 
the following form: 



S Rcs {x)= dtexp(-t)t b [M/N}(xt). (8) 



Similar to the Pade-analysis, within the Pade-Borel- 
Leroy resummation procedure, various final estimates of 
a resummed series can be obtained depending on the type 
of the Pade approximants chosen. In addition the fit pa- 
rameter b can be used for adjusting the resummation 
procedure to provide a self-consistent convergence of the 
results. In Fig. we present the result of the e-expansion 
processing by means of the resummation procedure ((?])- 
(||). Here, the estimates of m c on the basis of the higher 
three-, four- and five- loop e-expansion are depicted de- 
pending on the approximant type as well as on the fit 
parameter b. Among all the possibilities we chose the 
approximants close to the main diagonal, namely [1/1], 
[0/2] in three-, [2/1], [1/2] in four and [3/1], [2/2], [1/3] 
in the five-loop approximation. Notably, approximants 
having poles on the real positive semiaxis are considered 
as unreliable. For instance, Fig. || displays the real part 
of values with imaginary parts smaller than 10 -7 . 

Since no information is available on the large-order be- 
haviour of series (^) , all the values of b can be considered 
as equally suitable. An evident property of all the curves 
in Fig. ^| is a saturation of the value of m c for large val- 
ues of b. In order to exclude the values on the basis of 
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large b and taking into account that in similar tasks of 
critical phenomena b accounts for several units, we re- 
strict ourselves only to the values of b within the interval 
< b < 5. We neither take into account the approxi- 
mant [1/2] since it provides only two values within the 
interval. 



2.0 



1 .8 
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der of perturbation theory we take up the values on the 
basis of different approximants of the same order as in- 
dependent. Doing so allows us to consider their average 
as an overall estimate of m c within a given number of 
loops. The overall estimates of m c obtained in this way 
depend on the loop order oscillatively, which permits to 
suppose that a successive central value lies between two 
preceding central values. Thus we average the overall es- 
timate on the basis of the approximation order with a 
corresponding overall estimate from the preceding order, 
choosing the error bars as a half of a difference between 
maximal and minimal values from these pairs. Finally, 
we obtain a sequence of confidence intervals of m c such 
that error bars of a higher result lie completely within 
error bars of the previous one. 

Based on the analysis described above we obtain the 
following estimates of m c from the e-expansion (ph: 



4LA: m c = 1.996 ±0.104, 
5LA : m r = 1.923 ±0.051. 



(9) 



One can note from estimates (||) that only the five-loop 
e-expansion excludes the values m c > 2, and therefore 
brings about a weakly quenched disorder irrelevancy for 
the xy-model universality class. 



Fig. 2. The estimate of the marginal order parameter di- 
mension m c obtained on the basis of Pade-Borel-Leroy re- 
summation of the e-expansion for different values of the 
fit parameter b. The following symbols show the results based 
on different Pade approximants: ■ and • mean the [1/1], [0/2] 
approximants, ▲ and ▼ show [2/1] and [1/2] approximants, 
while ♦, + and x denote the [3/1], [2/2] and [1/3] approxi- 
mants respectively. 



To obtain a confidence interval for m c we analyse the 
data of Fig. || in the following manner. As a first step 
we obtain the central value of m c on falling look on the 
data on each approximant separately. To this end we av- 
erage all the values on the basis of each approximant 
over 6, < b < 5. The corresponding error bars are 
then determined by a half of the difference between the 
maximal and minimal values of m c within the consid- 
ered interval of b. The resulting estimates of m c in the 
successive number of loops are presented in Fig. [|, a, 
where the result on the basis of each working approx- 
imant is shown by its central value (black circles) and 
the error bars. For instance, in the three-loop approxi- 
mation two confidence intervals correspond to [1/1] (a 
smaller central value) and [0/2] (a larger central value) 
approximants; the four-loop approximation contributes 
by a single estimate on the basis of [2/1] approximant, 
while the five-loop approximation yields three confidence 
intervals from the [1/3] (the largest central value), [3/1] 
(the middle central value) and [2/2] approximants. 

In order to get a final estimate of m c within each or- 



1 .93 -, 




1 .88 I | i | i | i | i | i | i | i | i | i | i | i | i | i | i | i | i | i 

1 2 3 4 5 6 7 8 9 1 1 1 1 2 13 14 1 5 16 
b 

Fig. 3. The estimate of the marginal order parameter di- 
mension m c obtained on the basis of Pade-Borel-Leroy re- 
summation of the pseudo-e expansion ([f|) for different values 
of the fit parameter b. The following symbols show the results 
based on different Pade approximants: ■ A, T denote the re- 
sults on the base of [1/1], [2/1], [1/2] approximants, + and 
x show the [3/1], [1/3], [2/2] approximants, while * and — 
correspond to [2/3] and [3/2] approximants respectively. 

Let us turn now to the Pade-Borel-Leroy analysis of 
the pseudo-e expansion (^). The dependencies of the es- 
timates of m c on the basis of the pseudo-e expansion on 



MARGINAL DIMENSION OF A WEAKLY DILUTED QUENCHED m- VECTOR MODEL 



the fit parameter b are shown for different approximants 
types increasing number of loops from three to six in 
Fig. Since the pseudo-e expansion is one order longer, 
in addition to the approximants of Fig. ^ we use here 
the approximants of the six-loop order, namely the near- 
diagonal approximants [3/2] and [2/3]. The values of m c 
obtained on the basis of the approximant [0/2] do not fit 
Fig. H and thus are not presented, though they will be 
taken into account in the final calculations. 

Applying the procedure identical to the e-expansion 
we obtain the resulting estimates of m c depending on the 
type of the approximant as presented in Fig. ||, b. For in- 
stance, in the three-loop approximation two confidence 
intervals correspond to [1/1] (smaller central value) and 
[0/2] (larger central value) approximants; the four- loop 
approximation contributes by two almost identical esti- 
mates on the basis of [1/2] and [2/1] approximants, while 
the five-loop and the six-loop approximations yield three 
and two confidence intervals respectively. Subsequently, 



averaging the overall estimate on the basis of a given 
loop order with the corresponding overall estimate from 
a preceding order, one obtains a sequence of convergent 
confidence intervals for m c : 



4LA: m c = 1.976 ±0.117, 
5LA: m c = 1.904 ±0.013, 
6LA: m c = 1.912 ± 0.004. 



(10) 



The sequence of results (|T(j) obtained on the basis of the 
pseudo-e expansion (|J) shows evidently better conver- 
gence properties than the corresponding results @ ob- 
tained on the basis of the e-expansion (0). For instance, 
the error bar of the five-loop e-expansion estimate is four 
times larger than the corresponding estimate on the basis 
of pseudo-e. 



number of loops 



number of loops 



Fig. 4. The confidence intervals for the marginal order parameter dimension m c of a weakly diluted m-vector model as 
functions of the order of approximation. The left picture presents the result of the Pade-Borel-Leroy resummation of the 
e-expansion ([]), the right-hand picture corresponds to the analysis of the pseudo-e expansion (Q). See the text for the whole 
description. 



III. CONCLUSIONS 

In the theory of critical phenomena the values of crit- 
ical exponents and critical amplitudes ratios have at- 
tracted considerable interest of researchers. This study, 
however, has calculated another important quantity, 
which is responsible for a crossover phenomenon in the 
weakly diluted quenched m-vector model. Based on the 
Harris criterion as well as on the field-theoretical renor- 
malization group approach we determined the marginal 
order parameter dimension m c of the model. Before our 
work only one reliable numerical estimate of m c was ob- 



tained If]. The result m c = 1.942 ±0.026 |§ stemmed 
from the application of a conformal mapping resumma- 
tion technique to the massive six-loop RG functions at 
d = 3. We completed this study with two more results 
from two different RG schemes. Within the minimal sub- 
traction scheme m c is obtained as an e-expansion (||), 
while within the massive scheme m c is calculated as a 
pseudo-e expansion (Q). Both series allow to estimate 
a numerical value of m c after their analysis by Pade- 
Borel-Leroy resummation procedure, however the final 
estimate of the paper m c = 1.912 ± 0.004 is based on 
the pseudo-e expansion. This is explained not only by 
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a longer series calculated for the expansion, but also by 
its better convergent properties in comparison with the 
e-expansion. Such a situation has been already observed 
for a calculation of a marginal dimension of a cubic model 

a 

Though the obtained value of m c is very close to the 
integer value 2, our estimate evidently shows that within 
the error bars it is slightly smaller than 2. This result im- 



plies that universal properties of the xy-model are not af- 
fected by a weakly quenched disorder at criticality. This 
conclusion is confirmed by recent experimental studies 
on the critical behaviour of the superfluid He 4 in porous 
media ||. 

We thank Claude Bcrvillier for useful discussions and 
Ted Knoy for a thorough reading of the manuscript. 
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